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Discrete collocation method for 
Volterra type weakly singular integral 
equations with logarithmic kernels 


P. Mokhtary* 


Abstract 


An efficient discrete collocation method for solving Volterra type weakly 
singular integral equations with logarithmic kernels is investigated. One of 
features of these equations is that, in general the first derivative of solution 
behaves like as a logarithmic function, which is not continuous at the origin. 
In this paper, to make a compatible approximate solution with the exact 
ones, we introduce a new collocation approach, which applies the Miintz- 
logarithmic polynomials(Miintz polynomials with logarithmic terms) as basis 
functions. Moreover, since implementation of this technique leads to integrals 
with logarithmic singularities that are often difficult to solve numerically, we 
apply a suitable quadrature method that allows the exact evaluation of inte- 
grals of polynomials with logarithmic weights. To this end, we first remind 
the well-known Jacobi—Gauss quadrature and then extend it to integrals with 
logarithmic weights. Convergence analysis of the proposed scheme are pre- 
sented, and some numerical results are illustrated to demonstrate the effi- 
ciency and accuracy of the proposed method. 


Keywords: Discrete collocation method; Miintz-logarithmic polynomials; 
Quadrature method; Volterra type weakly singular integral equations with 
logarithmic kernels. 


1 Introduction 


In this paper, we develop an approximate approach to obtain the numerical 
solution of the following Volterra type weakly singular integral equation with 
logarithmic kernel 
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y(x) = g(a)+ f in (a —t)K(a,t)y(t)dt, ce Q=([0,1], O<t<a<1, (1) 
0 


where the continuous functions g(x) and K(a,t) are given and y(z) is the 
unknown solution. Such kinds of equations arise from solution of Dirichlet’s 
problem for the Laplace equation in two dimensions in terms of a single-layer 
logarithmic potential [24], solution of the reduced wave equation in two di- 
mensions using an integral equation with a kernel, which can be expressed 
as a Hankel function of order zero that this kernel, has also a logarithmic 
singularity [25], investigation of electrostatic and low frequency electromag- 
netic problems [15], methods of computing the conformal mapping of a given 
domain [28], solution of electromagnetic scattering problems( [2], [6]), deter- 
mination of propagation of acoustical and elastical waves( [5], [6]), boundary 
value problems of plane elasticity theory for regions with a defect [3], prob- 
lems of diffraction by thin screens [12] and so on. 


Weakly singular integral equations with logarithmic kernels are usually 
difficult to solve analytically; so it is necessary to provide reliable numerical 
techniques. There are several approximate methods proposed to obtain the 
numerical solution of these types of equations, which we refer to some of them. 
In [4], authors designed a computational meshless discrete Galerkin method 
to solve the second kind Fredholm integral equations with logarithmic kernels. 
In [17], authors developed a collocation method for the numerical solution 
of a special integro-differential equation with logarithmic kernel using airfoil 
polynomials of the first kind. A collocation method based on the periodic 
splines was introduced in [LI] to solve some logarithmic kernel integral equa- 
tions on open arcs. In [29], authors studied a special integral equation with 
logarithmic kernel and solved it using product integration method. In [14], 
a piecewise Chebyshev expansion was considered to solve Volterra integral 
equations with logarithmic singularities in their kernels. The properties of 
the integro-differential equations of the convolution on a finite interval with 
kernel having a logarithmic singularity were studied in [8]. In [1], authors 
investigated two numerical approaches by means of an analytical integration 
in the vicinity of the singular point and extraction of the singular part. A 
Gauss type quadrature method with a logarithmic weight function was ex- 
tended to evaluate of the Cauchy type integral equations with logarithmic 
kernels in [BU]. 


In this paper, we design and analyze a reliable discrete collocation tech- 
nique to obtain a suitable approximate solution of (). Our strategy is based 
on the following two stages. From well-known existence and uniqueness the- 
orems [8], we can see that the first derivative of the exact solution has a 
singularity at the origin and behaves like as a logarithmic function. Then 
to establish a highly accurate approximate solution, it is necessary to rep- 
resent the collocation solution as a linear combination of the suitable bases 
functions with logarithmic asymptotic behavior like as Miintz-logarithmic 
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polynomials [8]. Since in the implementation procedure, integrals with log- 
arithmic singularities are observed to highly accurate evaluation of them, we 
use a generalized Gauss type quadrature with logarithmic weight function 
that calculates exactly integrals of polynomials [Z]. 

The reminder of this paper is organized as follows. In the next section we 
present the required preliminaries for our subsequent development. Here, we 
introduce the Miintz-logarithmic polynomials as well as the generalized Gauss 
quadrature method for the integrals with logarithmic weight functions. In 
Section 3, we explain the application of the discrete collocation method using 
Miintz-logarithmic polynomials to approximate the solution of (M0). In Section 
4, we provide a reliable convergence analysis for the proposed algorithm that 
justifies the L?—norm of the error function tends to zero as the approximation 
degree tends to infinity. Section 5 devotes to our numerical illustrative and 
Section 6 contains our conclusions. 


2 Preliminaries 


In this section, we give some preliminaries that are required in the sequel. 


2.1 Miuntz-Logarithmic polynomials 


This subsection is devoted to a brief introduction on the Miintz-logarithmic 
polynomials. All of the details presented in this section a long with further 
details can be found in [[8j. 

The Mintz-logarithmic polynomials are defined as 


M,y(£) = Rn(x) + S,(x) nz, n>0, c EQ, (2) 
where R,,(x) and S,(x) are algebraic polynomials of degree [#] and [45+], 


respectively; that is, 


[24] 
Ry(x) = rv’, S,(x) = S- 3a". 
i=0 


It is shown that these polynomials are orthogonal with respect to the 
weight function w(x) = 1. Explicit expressions of the coefficients are obtained 
as follows: 


Theorem 1. 


e Ifn=2m, m= 0, we have 
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27 \2 an 
2m+1 . 2741 
( i Ea +2(m—i) 20 Gaara |) 
J=0,j Ai 


2 2 
pa (™PP) (™) | amet ye 
n= { m ) ("") es a aac oe G—-y)G+i+1) |? 
+4 . m : 
s=(m itt) (™ ‘) @ , O<i<m 
m i 


Proof. See [3]. 


The first few Miintz-logarithmic polynomials are given by 


Mo(a) 

My(2) 

Mo(a) = -3+ 42 -Inz, 

M3(x%) = 9 — 8a 4+ 2(14+ 6x) Ing, 
) 
) 
) 


Ma(x) = —11 — 24x + 36x? — 2(1 + 18x) Ing, 
Ms(ax) = 19 + 276 — 294x + 3(1 + 48a + 6027) Ina, 
Meo(x) = —21 — 768 + 390x? + 400x* — 3(1 + 96x + 300x”) Ina. 


2.2 Jacobi-Gauss quadrature 


This subsection presents the mechanism of the Jacobi-Gauss quadrature [[Z, 
9,27], which is to seek the best numerical approximation of an integral by 
selecting optimal nodes at which the integrand is evaluated. The Nth-order 
Jacobi-Gauss quadrature formula (7G)%" to approximate the integral is 
given by 
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t N 


i f(r) — r)(1 4 rr & (TE)"8(f) = Yo Flat we, 


Ae i=1 
a,8>-—-1, reé[-1,]], 


which the integral is evaluated exactly if f(r) is a polynomial of degree 2N — 
1 or less. The nodes {2° “31N | and the corresponding weights (ee 
depend on a, 3 and are given by the following formulas [1,27]: 


e {xo B }W_| are the zeros of the following orthonormal polynomial 


nN ()=\laraHT(N Fa+DI(N+ A841) ON 


where I (r) is the classical Jacobi polynomial of degree N. 


e The weights {W; BIN | are given by 


N-1 


2: 
“= (Peter ). 1<i<N. 


n=0 


The error term for the Nth-order Jacobi-Gauss quadrature (7G)~? is 
given by [Zj 


1 (2N) 


(3) 


-1 
where € lies somewhere on [—1, 1] and 


2N+o+6+1NIT(N +a+1)0(N+6+1)P(N+a+641) 
QN+a+B+1)TQN 404841)? 


ay? = 
2.3 Gauss type quadrature for logarithmic weights 
integrals 


This subsection devotes to introduce a generalized Jacobi-Gauss quadratures 
(Z] to approximate the following integrals 
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(a) / f(r)In(1— (1 — nar, 


—1 
1 


(b) io In(l+r)(1+ r)edr. 


—-1 


Considering the first integral in (@) we have 


1 1 

0 a _ 0 aln(1—r) 

a fan f(r)dr = a In(l—r) F(r)dr 
—1 -1 


1 


= [© In(l—r)(1—1r)dr, 


-1 
and equivalently 


1 


[fina -na-n"ar = je f ans (rar 


—1 


Dd 
~ 


a 


da 
= (GIG) (f), 


& da 


N 
(Ja) = (vee 


a,0 a,0 
(= fo) + wee 
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yu) 
ree") 


where “GG” is an abbreviation for the generalized Jacobi-Gauss. From 
(B), we can obtain the following formula for the error term of Nth order 


generalized Jacobi-Gauss quadrature (G IG)~: 
1 
(fina are ner) — (GIG)*"(f) 
=i) 


= F Ewalo, 0, f(r)) = Enr(a, 0, f(r)) 


1 
| WS 
(909 IN+atl 


+2U(N+a+1)-2UQ2QN+a+ D )exla 0,f(r)), (6) 
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where U(z) = 2” Me) — Te is the psi or digamma function. 


Proceeding the same technique with ()-(@) we can obtain the following 
approximation for the second integral of (@) 


1 


: 
n r tps r)? f( 
fro (l+r)(1+r)’d a5 | (l+r)’ f(r 


-1 


dW; 0,8 “soda 1/,0,8 
Eebmn )+W; f(z; ae 


= (GIG) (f), 


which has the following error function 


1 


(f rea +r? in tnar) - Eze) 


-1 


= spene(O B, f(r)) = En,(0, 2; f(r)) 
={In(2 
= Me) Sir ai 


+2U(N +68 +1)—2WQN + 8+ Denso B,f(r)). (7) 


The relations (@) and (@) conclude that the generalized Jacobi-Gauss 
quadratures (GIG)°(f) and (@7G)°%(f) calculate the integrals of (@) 
exactly for polynomials of degree 2N — 1 or less same as the Jacobi-Gauss 
quadrature (7G)“9(f). 


3 Numerical approach 


The main concern of this section is to obtain the discrete collocation solu- 
tion of (M) when the Miintz-logarithmic polynomials are applied as the basis 
functions. To this end we represent the collocation solution of (I) as 


N 
= »; anMn(2), (8) 
n=0 
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such that the unknowns {a,,}+_, satisfy in the following linear algebraic 
system: 


vi 


in aaa jo (e; — t)K(ai,t)yn(t)dt, O<i<N, (9) 
0 


where {x;}/9 is the shifted Legendre-Gauss quadrature points [27]. Apply- 
ing the variable transformation 


t=+,;(0) = 5 (0+ 1), 6 ¢€ [-1,1], 


we can rewrite (Q) as follows: 


1 
+ fima- ko Dunieoyao} (10) 
=i 
for 0<i<N and K(a;,0) = K(a;,t;(0)). Substituting (8) into (0) yields 


N 
> oof Mate - 4 nay) rage =g(ai), 0<i<N, (11) 


com / In (1 — 0) K(2r;,0)M(t(6)) a9. 


Using (2), we have 


Discrete collocation method for Volterra type weakly singular ... 103 


Afr) = [Ke8) (Rata) + S,(t;()) In (n() ) a0 


where 
FL) (0) = K (ai, 0)Sn(ti(6)). (12) 


Using quadratures (7G)%* and (G7G)°-", we obtain 


Ayn? ce AM? = (TG)O(FL (0) + | (GTE)* (Fy"*? (6) 


On the other hand, from (Q) we can write 


AG) = | In (1 — 0) K(a;, 0) (R (t;(0)) + S,(t;(0)) In (u(0))) do 


1 
In (1 — 0) FA") (0)d0 + i In (1 — 8) In (1 + 0) FS" (0) d0 
—-1 


lI 
= en pei 


1,n,2x;) Ale) (14) 


where 
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AS) = | inti — 0) Fi" (9) a0, 


AQ) — | In (1 — 8) In (1 + OFS"? (8)d0. 


| 
| 


Trivially we can write 


ln, x; 1,n,x% 
A, laid Ag N= 


~ 


(GIG) FY™ EO), (15) 


and 


0) 
Agr?) = fin —e)in( + VF (as 
Si 
1 
% / In (1 — 6) In (1 + 0) FS" (0) 0. (16) 
0 


Using the variable transformation 6 = 0)(s) = 
and @ = —6;(s) in the first integral of (I8), we obtain 


s+ in the second integral 


re 1 l—s 34+8 NL; NL; 
Apne) a1 Piglet S#8 (ze oy(ay + 70% (-000)) a 


-1 


t 1 

iL N,Xy Ni,XL; 

5 [mao , Y(s)ds—In2 fA ou, 
Zi 


-1 


where 


Fee) =n 2 F2 (Amal +7 —Hi9). AN 
then we can write 


; ; 1 
2.n,0i) | 2,n,24) __ 
Apres AN = 5 

a=0 


([ezarreren] 
- maze) . (18) 


Substituting (8) and (@) into (G) yields 
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TBs NL; 1n,x; 2,1, 2X4 
AP OSA A PA (19) 


Inserting (9) and (23) into (ZU), we can conclude that the discrete collo- 
cation solution of (M) is characterized by 


where the unknowns {@,,}/_9 satisfy in the following system of linear alge- 
braic equations 


The matrix formulation of (20) is given by 


aM = g, (21) 


where a= [G@o, @1,...,4y], C= [9(@o), 9(#1),-- .,g(an)]*, and 


M= My, = (AiDi + A2D2), 


such that 
Mo(x%o) Mo(21) Mo(xw) 
re My,(%o) Mi(21) «+> Mi(zy) 
l= . . . : ; 
My (xo) Mwy (a1) -+- Mn (an) 
si) a an 
, Ain’ AN’ ae Ain” 
1— F é . ’ 
Afheo Age? agen 
ADN? Ag? Ag’? 
1,2 1.x l,a 
ae An A as A 
2 — , F é ’ 


AG A ar Ae 


and D,; and D2 are the diagonal matrices with diagonal entries 
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4 Convergence analysis 


In this section, we provide a reliable error analysis for the proposed technique 
to justify convergence of the proposed approach. In our analysis we shall ap- 
ply the following definitions and lemmas. 


Definition 1.[ [I], (9-23), [27] 
1. £0) = {u) ful = flu(a) Pa < oo} 


2. C(Q) is the space of all continuous functions on (. 


3. Iyu is the Legendre—Gauss interpolation polynomial and defines by 


N 
Inu(x) = 5 u(a;)Li (x), 
i=0 
where L;(x), 0 <i < N, are the Lagrange interpolation basis functions 
associated with the Legendre—Gauss points {x;}_o. 


4. Let (4, ||.||) and (J, ||.||) be normed vector spaces, and let K: ¥ > Y 
be a linear operator. Then K is compact, if the set 


{Kul |lulla < 1} 


has compact closure in Y. For example, the integral operators with 


continuous and weakly singular kernels are compact operators on C(Q) 
and L?(Q). 


Lemma 1.[see [R]] Let K(x,t) € C(Q x Q); then for any g(a) € C(Q), the 
Volterra type weakly singular integral equation with logarithmic kernel ()) 
possesses a unique solution y(a) € C(Q). 


Lemma 2.|see [I] Let ¥ be a Banach space, and let KK: X + X be compact. 
Then the equation (A—K)u= f, AZ#0 has a unique solution u € X if and 
only if the homogeneous equations (\—K)z = 0 have only the trivial solution 
z = 0. In this case the operator \—K : ¥ > &X has a bounded inverse 
(A—K)71. 


Lemma 3.|see [26]] [f u(x) € C(Q), then we have 


lu —Inull2 70, as N- ov. 
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Lemma 4.|see [I9]] For every bounded function u(x), there exists a constant 
C independent of u(x) such that 


sup ||Ivull2 < Cllulloc- 
N 


Theorem 2. Under assumptions of Lemma 1, assume that the following 
conditions are satisfied: 


1. The given functions g(x), K(a,t), and the approximate solution 9x (x) 
are continuous on their domains. 


2. The homogeneous equation (Dl) with g(x) = 0 has the only trivial solu- 
tion. 


3. §n(x) € C(Q) and the numerical integral operator 
N : 
castn Sone {nah +a} 


is a bounded operator on C(Q) to C(Q) with Ayn and Aj’y, defined in 
(13) and (19), respectively. 


4. The quadrature errors 


Ey,0(0,0, Fy" (9)) loo: llEw,s (0,0, F3°" ()) loos 


[éw,0(0, 8, F3°*(8))| 


[éw.0(@, 0, F7*(6))] sally 


= 
coll 


léw.s(0, 0,F3'*(s))| O<n<QN, 


9 
wall. 


converge to zero as N + co. Here the functions Fy" (0), Fy" (8) and 
F3°"(s) are given in (TQ) and (Mi), respectively, and the error terms 
En,r(0, 0, FW), a 

En,r(a, 0, f(r)), En.r(0, 8, f(r)) are defined in (B), (G), and (@), respec- 
tively. 


Then we have 
Nim |ly(x) — gn (@)|l2 = 0. 
—0o 


Proof. Using (20), we can conclude that the discrete collocation solution 
Yn (a) for the equation (f) satisfies in the following operator system: 


Yn (#i) = g(@i) + Kn yn (ai), OStSN, (22) 
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where the operator Ky is the numerical integral operator defined in the 
assumption 3. Multiplying L;(x) on both sides of (22) and summing up from 
i=0toi=N yield 


In (9n) = Ing + InKngn.- (23) 
Subtracting (M) from (23) gives 
y — In(9n) = (9 — Ing) + Ky — InKngy), 
or equivalently 
(id — K)én = ery (g) — ery (Yn) + ery (K9N) + In (Kin = Kyi), (24) 


where €y = y(x) — Yn (2) is the discrete collocation error function, e7,,u = 
u — Inu is the interpolation error function, id is the identity operator, and 
K is the following continuous and compact integral operator with weakly 
singular logarithmic kernel 


Applying Lemma 2 with ¥ = C(Q) and using the assumption 2, the 
relation (24) can be rewritten as 


llénll2 < id — KJ" oo (er (g)ll2 + llery Gn )Il2 


+ |lery (KY lo + IL (Kn — Kuan) (25) 


with ||(id — K)~*||,. < oo. Due to Lemma 3 and assumption 1 we have 
llery(g)ll2, llery Gn lla, llery(KYw)|2 40, as Noo. (26) 


Now, it is sufficient that we show ||Iy (Kyn — Kngn) le > 0as N > 00. 
To this end, using Lemma 4 and assumption 3, we can write 


I| Iw (Kon — Kn Yn) |l2 < Callow — Kw Gnllo; (27) 


which Kyy — Kyyy is the numerical integration error function. According 
to the definition of Ky) and numerical approach proposed in the previous 
section, we can deduce that the numerical integration error is established by 
calculating the errors obtained from the applying Jacobi-Gauss and general- 
ized Jacobi-Gauss formulas in (3), (3), (8), and (1). Consequently, we 
can write (27) as follows: 
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In (Kon — Kyi) |lo < Co 3 (ew 0, Fy" (8))Ihoo 
n=0 
-[[éu0a79"0), 
a: [évo(a,0.F*(0))] | 
a [évs(o,0.F3%(8))] 


+ Exl0. 0.79"). 


and thereby using assumption 4 in the inequality above, we deduce 


In (Kyn —Knyn)|l2 70, as NO. (28) 


Finally, the required result can be obtained by applying (26) and (28) in 
(25). 


5 Numerical Results 


In this section, we illustrate some examples using the method proposed in the 
previous section and confirm its validity. All of calculations performed on a 
PC running Mathematica software. In the obtained results we presented some 
essential items regarding the L?-norms of the error functions and compari- 
son results between our scheme and the well known Chebyshev collocation 
method [9,27]. 


Example 1. Consider the following problem 


x 


yx) = g(x) + f In (e—t)y(tat 
0 
with ‘ 
g(x) = z(Inxz —1)+ oo (7 —21+18In2 —6l1n? c) 
and y(#) = x(Ina# — 1) as the exact solution. 


We solve this problem by the proposed method with N = 3. Indeed, we 
seek a discrete collocation solution in the form 
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to approximate this problem and find its unknown coefficients such that they 
satisfy in the linear system (21) with N = 3. Considering the Gauss—Legendre 
collocation points 


xo = 0.330009, x; = 0.669991, x2 = 0.0694318, x3 = 0.930568, 


and the operational matrices 


1 1 1 1 
MM, — { 70-1086 0.5995. —1.6674 0.9280 
1= | 9.5713 0.0805 —0.0549 0.7942 | ° 
0.2477 —0.3808 0.8873 0.6080 
2 2 
ra re nets 5.3348 —0.1439 
1= | 0.4627 0.5191 1.6126 —0.1338 | ° 
0.5550 —0.2017 0.1359 —0.1194 
0.6137 —0.6137 —0.6137 —0.6137 
a, — | 79-6095 -1.0441 0.3471 -1.2457 
2=]| 0.7718 0.1092 0.6562 —0.5302 |’ 
0.1912 0.6068 —1.3532 0.0376 
0. 2073 0 0 
a: se 0 0 
ae —0.1167 0 ; 
0 —0.3560 
0.1650 0 0 
0 0 
a te o |? 
0 0.4653 


the system of linear algebraic equations (21) is presented in the following 
form 
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1.0449 + 1.6959a9 — 0.66734, — 0.836242 — 0.1142a3 = 0, 
1.6603 + 1.9383a9 + 0.65584, — 0.146342 — 0.6580a3 = 0, 
0.2955 + 1.254649 — 2.30194, + 0.1105a2 + 0.9501a3 = 0, 
1.8965 + 1.9975a9 + 1.4564a, + 0.993342 + 0.5830a3 = 0, 


which its solution gives by 


Gp = —0.75, a =—0.25, do = —0.0833, a3 = 0.0833. 


Consequently the discrete collocation solution #3(a) represents by 


93(%) = —5.2953 x 10-7 — w — 1.3027 x 107” na + 0.99992 Ina. 


The L?-norm of the error function is 1.3445 x 107!° that is in a very 
good agreement with the exact ones whereas the approximation degree is 
very small(N = 3). 


Example 2. Consider the following problem; 


x 


y(a) = g(a) + / In (a — te* ty (t)adt, 
0 


with ; 
g(z) =e" *Inzet+ =(- 124+ 77 — ina(—2+In2)). 


and y(z) = e~* Ina, as the exact solution. 


We solve the problem and report the obtained results in Table [] and Fig- 
ure I. In Figure I], we plot the L?-norm of the error function in terms of 
the various values of the degree of approximation N. Figure [I] shows that 
the proposed algorithm obtains a good accuracy with suitable values of NV. 
Moreover to make a comparison, we also solve this problem by implementing 
the well-known Chebyshev collocation method [9,27] and give the obtained 
results in Table I. Based on Table 0], we confirm that the Muntz-logarithmic 
polynomials makes faster rate of convergence for the discrete collocation so- 
lution of this problem compared with the classical Chebyshev polynomials. 


Example 3. Consider the following problem: 


y(a) = g(a) +2 / In (« — #) ty (t)dt, 
0 
with m 
222 (—13016 + 6930 In 2 + 3465 ln x) 
38115 : 
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Table 1: The numerical errors of Example 2 


Numerical Errors 


N Our Method Chebyshev collocation Method [9,27] 
2 3.61 x 107° 1.3 x 107! 
4 8.13 x 1076 4.22 x 107? 
6 9.38 x 10-8 2.14% 10-7 
8 9.84 x 10719 12% 10-* 
10 3.44 x 1071? 8.73 x 1073 
12 6.37 x 10718 6.28 x 1073 


Log,,( Error) 


Figure 1: Obtained L?— norm of the error function versus N for Example 2 


and y(x) = x7,/z as the exact solution. 

We have calculated the approximate solution with different values of N 
and displayed the obtained results in Table 2 and Figure 2. Table 2 and 
Figure Q, present the L?-norm of the error functions versus N. As it can be 
observed although the exact solution has singularity at zero, the numerical 
errors are decreased with an appropriate rate as the approximation degree N 
is increased. 

Example 4. Consider the following problem: 


x 


y(x) = g(x) + / In (a — t)y(t)at 


0 
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Table 2: The numerical errors of Example 3 


N Numerical Errors 
2 2.69 x 10-2 
4 3.32 x 10-4 
6 1.18 x 1077 
8 1.4 x 10710 
10 1.84 x 10-7 
12 5.43 x 10715 


Log,,( Error) 


Figure 2: Obtained L?— norm of the error function versus N for Example 3 


with 
g(x) = e*(1 + y+ T(0,x)) +lInz, 


where 7 is Euler’s constant with the numerical valuex 0.577216 and ['(0, x) 
is incomplete gamma function satisfies 


co 


T(a,z) = pectetae 
z 
Here the exact solution is given by y(x) = e”. 


The obtained numerical results from implementation of the proposed dis- 
crete collocation scheme are reported in Table Bl and Figure B. The presented 


114 P. Mokhtary 


Table 3: The numerical errors of Example 4 


N Numerical Errors 
2 6.08 x 10-3 
4 2.01 x 10-4 
6 2.68 x 10-6 
8 1.78 x 1078 
10 6.99 x 10714 
12 1.85 x 10718 


results confirm that our scheme provides reliable results for smooth solutions 


of (O). 


Log,,( Error) 


Figure 3: Obtained L?— norm of the error function versus N for Example 4 
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6 Conclusion 


In this article, we developed a new discrete collocation method based on 
the Miintz-logarithmic polynomials as basis functions. Moreover, we used 
highly accurate Jacobi-Gauss and generalized Jacobi-Gauss quadratures to 
approximate the integrals with Jacobi and logarithmic weights, respectively. 
Convergence analysis of the proposed method were presented and some nu- 
merical examples were illustrated to confirm the applicability of the presented 
discrete collocation scheme. 
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